Distribution of level curvatures for the Anderson model at the 
localization-delocalization transition 
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We compute the distribution function of single-level curvatures, P{k), for a tight binding model with 
site disorder, on a cubic lattice. In metals P{k) is very close to the predictions of the random-matrix 
theory (RMT). In insulators P{k) has a logarithmically-normal form. At the Anderson localization- 
delocalization transition P{k) fits very well the proposed novel distribution P{k) oc with 
« 1.58, which approaches the RMT result for large k and is non-analytical at small k. We ascribe 
such a non-analiticity to the spatial multifractality of the critical wave functions. 

PACS numbers: 71.30.-)-h, 72.15.Rn, 05.60.+W 



An important consequence of the one-parameter scal- 
ing theory 0] of the Anderson transition is the existence 
of three universality classes for the energy level statis- 
tics of disordered systems. In the metallic regime they 
are described by the strongly correlated Wigner-Dyson 
(WD) distribution, while in the insulating regime they 
follow an uncorrelated Poissonic law. This difference has 
its fundamental origin in the underlying nature of the 
corresponding eigenstatcs, being extended and strongly 
overlapping in the first case but localized in the second. 

The properties of the third universal statistics, de- 
scribing the spectral correlations at the critical point, 
have been only recently the subject of intense investi- 
gation, both analytical and numerical. Using the re- 
sults of numerical simulations, Shklovskii et al. sug- 
gested that the spacing distribution function P{s) has 
the WD form P{s) ~ s for small s and the Poissonic 
tail P{s)s^ for large s. Further analytical investiga- 
tions showed that the two-level correlation function 
R{s) in the critical region has a novel power-law asymp- 
totic decay R{s) = —c/s'^~'^ with a nontrivial exponent 
7 = 1 — 1/ vd. Here v is the critical exponent of the cor- 
relation/localization length ^, which depends on the di- 
mensionality d. Thus the two-level correlation function in 
the critical region resembles qualitatively the WD func- 
tion which applies to the metallic phase. On the other 
hand, the level number variance Tj2{N) = {{SN)'^) in an 
energy strip of width A^A, (A is the mean level spacing), 
still contains a dominant Poissonic term, linear in 

N, which is typical for insulators. 

It can be shown that the Poissonic term in S2(^) 
is only possible if a normalization sum rule on R{s) is 
violated in the thermodynamic limit (TL) L — > c», (L is 
the system size) which signals a qualitative change in the 
statistics of wave functions ^'(r). An analogous situation 
occurs in certain random matrix ensembles with broken 
unitary symmetry, that turn out to describe the critical 
statistics E2(7V) and P(s) very well, g 

Since the Poisson statistics describe localized states, 
it was put forward [B[ that the cause of the linear term 



in 'E2iN) is the existence near the Anderson transition 
of pre- localized states, for which sharp peaks in |5'(r)p 
contribute significantly to the normalization integral. In 
fact the existence of the pre-localized states - even in 
good metals - had been already conjectured in Ref. px| ] 
for the interpretation of a slow current relaxation in dis- 
ordered conductors [Q. Recently this conjecture has 
been confirmed in a more direct way using the super- 
symmetric sigma-model. [|ll| , ^ . Furthermore, the mul- 
tifractal structure of wave functions ^'(r), expected at 
the Anderson transition, ||l^ can be represented as 
a superposition of peaks in |'I'(r)p corresponding to pre- 
localized states with a certain distribution of amplitudes 
and exponents for their power-law spatial decay. 

A statistical property that may shed some light on 
the problem of pre-localized states is the distribution of 
curvatures of the single energy- levels. The level curva- 
ture measures the sensitivity of the energy spectrum to 
a change of the boundary conditions. For each single 
energy level e„((/5), the curvature Kn is defined as 
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The parameter if enters in the definition of general- 
ized boundary conditions in the z direction satisfied 
by the single-level wave functions, ^'(a;,y, z + L) = 
e'''^^'^'^{x,y, z). For a quasi-onedimensional sample, the 
parameter if is equivalent to the phase generated by an 
Aharonov-Bohm flux 4> = ip(j)o that pierces a closed ring. 
The way in which an eigenvalue responds to a small twist 
of the boundary condition is obviously related to the 
nature of the corresponding wave function. Typically 
an eigenstate extended through the sample will feel any 
change in the boundary conditions and will have a large 
curvature, whereas a localized state, far from the edges, 
will be insensitive to the twist and its curvature will be 
close to zero. However we will see that this intuitive ar- 
gument must be taken with caution. 

In this paper we investigate numerically the curvature 



1 



distribution of the 3d Anderson model. Similar calcula- 
tions have been performed earlier by Zyczkowski et. al, 
psf , who studied the distribution for all the three regimes 
of the disorder, and by Braun and Montambaux who 
just looked at the metallic (diffusive) regime. The pur- 
pose of our work is to carry out an accurate scale anal- 
ysis of the full curvature distribution right at the metal- 
insulator transition, which has not been done so far. We 
will first compute very accurately the curvature distribu- 
tion for the metallic and insulating regime, and provide 
a physical interpretation of the numerical results in these 
two cases, based on the nature of their wave functions. 
We will then show that the critical distribution can be 
fitted extremely well by a functional expression similar 
to the RMT result that applies to metals, except for the 
presence of a new nontrivial exponent, which makes the 
distribution nonanalytical at A' = 0. We will finally ar- 
gue that such a behavior support the hypothesis of pre- 
localized states in the critical region. 

We consider a tight binding model on a cubic lattice 
of volume V = (L)^ , where L is is the number of sites in 
any direction. The one-particle Hamiltonian is 



H 



Y,e4c,+tYy' 
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The site energies Ci are randomly distributed with uni- 
form probability between — W/2 and -\-W/2. The param- 
eter W controls the amount of disorder in the system. 
The phase shifts Oij in the hopping term are chosen in 
such a way that when a particle hops from z to z + L it 
picks up a total phase shift equal to e*^'^'^ coming from 
generalized boundary conditions in the z direction. The 
Anderson transition is known to take place in such a sys- 
tem for Wc/t w 16.5 when (p = 0. 

For each disorder configuration we can compute the 
zero-flux curvature of each single level. Since in the 
presence of the disorder the level velocity u„ = 1 
is zero, the curvature can be calculated numerically in 
terms of the following finite difference 



Kni^ - 0) = 2 



[£„(^)-e„(0)] 
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In order to get accurate approximations of the differen- 
tials with finite differences, particular care must be taken 



in the choice of the ip. |15 1^ We fitted the energy shift 
[e„(±(/3) — e„(0)] to a parabola Kn^p"^ and checked that 
the value found for K„ is the same when we double ip. 
The optimal value of ip ranges form 10~^ to 10~^ de- 
pending on disorder and system size, increasing with dis- 
order. Our ensemble averages typically include several 
thousands disorder realizations, each one retaining half 
of the spectrum in an energy- window centered at e = 0.. 

We start with the metallic case where we know that 
classical RMT is applicable. Delande and Zakrewski p^ ] 
conjectured that the curvature distribution for complex 
systems described by RMT follows the expression 



(1 + A:2)l+/3/2 



(4) 



with P = 1,2,4 for the orthogonal, unitary and symplec- 
tic ensemble respectively. The constant Cp is fixed by 
normalization. The dimensionless curvature k is given 
by fc = K/ (K), where {K) is the first moment of the dis- 
tribution. It has been shown quite recently that 
Eq. is in fact exact for ensembles of large Gaussian 
matrices for all three symmetries. 
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FIG. 1. Distribution of curvatures in the metallic regime 
Wjt = 5, for different L, plotted with Eq. (|) for /3 = 1 as a 
comparison. In the inset P(log |fc|) = |fc|P(|fc|) vs log |fe|. 

In Fig. we plot our results for the distribution of 
the rescaled curvatures P(k) for several system sizes and 
W — ht < Wc- It is seen that the RMT expression, 
Eq. (^) with /3 = 1, is remarkably accurate in the metal- 
lic regime for the whole range of k. It is universal, in 
the sense that, after rescaling K, P{k) is scale invari- 
ant and disorder independent. Our results agree with 
Rcfs. p5| , p^ and settle the question of the 1/fc^ depen- 
dence of P(k) at large k, implying the divergence of the 
second moment. We emphasize that, for the 3d Ander- 
son model in the metallic regime, P{k) agrees very well 
with Vi3{k) of Eq. also near fc = 0, in contrast to 
the case of other quantum chaotic systems that dis- 
play nonuniversal features at small k, due to scarring of 
their wave functions. This point will be very useful in 
the interpretation of the critical P{k) near fc = 0. 

An important property of 'Pp{k) is the maximum at 
fc = 0. Intuitively, on the basis of the above argu- 
ments, one would expect that in the metallic phase, 
where the eigenf unctions are extended, most of the cur- 
vatures should be large. In fact, as we show below, the 
maximum in P(fc) is related to the "ergodic" nature or 
lack of specific internal structure of the typical metallic 
eigenstate. Using second-order perturbation theory, the 
second derivative of the eigenvalues with respect to the 
parameter if due to a small perturbation is written as 
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|(m|J^|n)| 
e™(0)-e„(0) 
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where |m) and 6^(0) are the eigenvectors and eigen- 
values of H at zero flux, respectively. is the cur- 
rent operator in the z direction. If we assume an or- 
dered spectrum — Cn = A(m — n), then for each 
term of this sum relative to the eigenvalue em+, sep- 
arated from e„ by an energy shift +5e„ = e„i^ — e„, 
there will be another term corresponding to em_ such 
that em_ — = —Sen- Moreover the matrix elements 
of these two terms |(m+| J2|n)| and |(m_|Jz|n)| must be 
equal if the corresponding unperturbed eigenstates are 
extended and structureless. Therefore we conclude that 
there is a cancellation in the sum of Eq. (^) , causing the 
curvature to be zero. Of course the actual distribution 
is not simply S{k) because the spectrum is not fully or- 
dered. Large values of the curvature |fc| ~ 1/s are due to 
anomalously small spacings s = |em(0) — e„(0)|/A <C 1 
for m = n + 1 or m = n — 1, which occur with probability 
P(s) ~ s^. Thus for large \k\ one obtains a power-law 
decay Vpik) J ds s>^ S{\k\ - 1/s) - |fc|-(2+/5) instead of 
zero, but the maximum of Vf}{k) remains at k = 0. 
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FIG. 2. Distribution of curvatures P(log|fc|) vs log for 
the insulating regime W/t = 30 for diflerent L's. The dotted, 
dashed, etc. lines are gaussian fittings. The solid line is the 
RMT result for P — 1 in these coordinates. 

This simple argument allows us to interpret also the 
results for the other extreme case, namely the regime 
of strong disorder where localization is present. Indeed, 
from Eq. (|) we see that now the cancellation in the sum 
does not occur: two states |m+) and I'm,-), separated in 
energy from \n) by ±(5e will have matrix elments that 
differ greatly from each other, since they are, in general, 
localized at different points. Therefore for a large but fi- 
nite system size L, the typical level curvature is small but 
non-zero. Thus P{k — > 0) = for an insulator. We can 
guess the form of the distribution function in this case if 
we assume that the curvature k is related to the ampli- 
tude of the typical wave function at the edges, log|fc| ^ 
log with l^fj ~ exp(— L/Lo), where Lo is the radius of 
the localized wave function. If the distribution of 1/io is 

a Gaussian P(l/Lo) oc e"*-""^-* centered around the in- 
verse localization length we come immediately to the 
conclusion that the curvature distribution for an insula- 



tor must be log-normal, P{\k\) = Cexp[— A(log |fc|-|-i3)^], 
with A ~ and B L/^. The numerical results of 

the simulations for the insulating regime are shown in 
Fig. |, where we plot P(log|fc|) = \k\P{\k\) vs. log|fc| 
for W = 30t > Wc and different system sizes. In these 
coordinates the distribution can be fitted rather well by 
gaussians, when the system size is large (L — 12), im- 
plying that P{k) is log-normal, in agreement with what 
is found in Ref. The vanishing of P(k) at fc = 

occurs very close to the origin and therefore is not easily 
resolved in a P{k) vs k plot. However for smaller L < 8, 
P(log |fc|) clearly deviates from a gaussian distribution at 
small k and approaches the RMT function even for such 
a large disorder. 
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FIG. 3. Distribution of curvatures P(\k\) vs |fc| at the criti- 
cal point {Wc/t = 16.5) for different L's. In the inset P{log\k\) 
vs log |fc| for L = 10, 12. 

We finally come to the central part of our work, the 
curvature distribution in the critical region. It was 
pointed out before that P{k) can be used to locate ef- 
ficiently the Anderson transition point, |jl^ since accord- 
ing to the one-parameter scaling P{k) must be scale in- 
variant at this point. However, there has been so far no 
detailed study of the full correlation function at the An- 
derson transition and its dependence on the sample size 
L. The results for different system sizes at Wc/t = 16.5 
are shown in Fig. 3 together with the RMT result V0=i (fc) 
for comparison. One can see that the curvature distribu- 
tion is close to Eq. ^ valid for metals. In particular, 
for large fc the two curves are indistinguishable within 
our numerical precision. Thus even at the critical point 
the P(fc) follows the |fc|~^ rule, which corresponds to a 
spacing distribution P{s) ^ s. On the other hand the 
curvature distribution P{k) at the critical point deviates 
slightly but significantly from Vp{k) at small fc. Here 
the maximum at fc = moves above the metallic limit 
and apparently is size dependent. This seems to indicate 
that, in contrast to the prediction of one-parameter scal- 
ing, even at the transition point there is a finite length 
scale R such as for L « R the full distribution is in- 
distinguishable from the RMT result, and for L >> i? it 
saturates to a complete scale-invariant curve that differs 
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from Eq. |^ at small k < l.As shown in Fig ^, where we 
plot the critical P{k) for L = 12, this curve fits very well 
the distribution function 



A 



(6) 



where A = /ir(3/^)[r(l/^)r(2/^)]-i > 1 is determined 
by normalization. Once normalized, the distribution in 
Eq. ^ obeys the condition (fc) — 1 for all iJ. = 2 — a. Thus 
the full distribution J-{k) is determined uniquely by its 
value at fc = 0, and with the known P{k = 0) the fit of 
the whole curve is parameter free. This procedure yields 
a = 0.42 ± 0.01. For further evidence we have plotted, 
in the inset of Fig. ^, the derivative of P{k) with respect 
to k, together with the derivative of ^(fc). Equation (||) 
implies that dP/d\k\ ^ for small 




P(|fc|) vs. |fc| at the critical point for L = 12. The 



dashed line is the fitting function of Eq. (f ) with a = 0.42. 
The solid line is the RMT function, Eq. (|), f3 = 1. In the 
inset dP(|fc|)/d|fc| calculated numerically from the P{\k\) data 
for L = 8, together with the derivatives of J-{k) and Vjiik). 

Such a non-analiticity in P{k) is the signature of the 
complexity of the wave functions at the critical point. 
Another (related) feature of the critical distribution P{k) 
is an enhanced probability to find a level with a small 
value of curvature. It can result from the multifractal 
spatial structure of the wave functions in the critical re- 
gion, provided that there is no drastic difference in ma- 
trix elements {m+\Jz\n) and (m_| J^jri) so that P{K — 0) 
remains nonzero. Indeed, such a structure can be con- 
sidered as a limiting case of "scarring" with a certain 
distribution of sharp peaks in |^(r)|^, which exhibits it- 
self in the spectrum of the multifractal dimensions dp < d 
seen in the inverse participation ratios (IPR) Ip- 



(7) 



Notice that in a metal Ip ^ \/L'^p~^^'^ and in the insula- 
tor Ip ^ constant in the TL. Thus one can say that the 
behavior at the critical point resembles the one in met- 
als, since in both cases the IPR vanishes when L ^ oo. 



However the exponents are different. This situation is 
analogous to what we find for the curvature distribution. 

In conclusion, on the basis of numerical diagonalization 
of the tight-binding Anderson model we have suggested 
a novel level curvature distribution function P(fc) at the 
Anderson transition, which exhibits a non-analytical be- 
havior for small k and approaches the RMT result for 
large k. We relate this non-analyticity to the spatial mul- 
tifractality of critical wave functions. 

One of us (C. M. C.) would like to thank R. Hlubina 
for interesting discussions. 
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